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O ■ ABSTRACT 

o ■ 

Solid particles with the composition of interstellar dust and power-law size distri- 
bution dn/da oc a"^ for a < Omax with Cmax ^ 3A and 3 < p < 4 will have submm 
^ . opacity spectral index /?(A) = dlnK/dlnu (p — 3)/?ism, where Pism ~ 1-7 is the opacity 

spectral index of interstellar dust material in the Rayleigh limit. For the power-law 
index p ~ 3.5 that characterizes interstellar dust, and that appears likely for particles 
\ growing by agglomeration in protoplanetary disks, grain growth to sizes a> 3 mm will 

Q>^ I result in /3(lmm) < 1. Grain growth can naturally account for /? ~ 1 observed for 

' protoplanetary disks, provided that amax ^ 3A. 

!>■ 

O _ 

■ Subject headings: planetary systems: protoplanetary disks dust, extinction, submillime- 

ter 
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|H !• Introduction 



From mm to mid-infrared wavelengths, the radiation emitted by protoplanetary disks is pri- 
marily to due to thermal emission from solid particles, just as for the interstellar clouds out of 
which the protostars form. If the emitting medium is optically thin, then the observed flux density 



a [ Fu is proportional to the mass of emitting particles, 

^ K{u)MMTd)d-^ , (1) 

where k(z^) and are the opacity and temperature of the solid material, Bi, is the blackbody 
intensity, and d is the distance. At submm wavelengths, the Rayleigh-Jeans limit hv <^ kT^ usually 
applies, and 

F^ « ^^{^)^p- (2) 

is proportional to the product M^Td- If can be estimated from the observed spectrum or other 
considerations, the observed F^, allows to be determined, provided that k and d are known. 

If the opacity has a power-law dependence on frequency, k(z^) cx z^^, the observed flux Fy cc v"^, 
with a = 2 + (3. Even if the absolute opacity cannot be determined, multiwavelength observations 
allow the opacity spectral index /3(A) = dh\ n/ dlxiv to be measured. 
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Observations of the infrared and submm emission from diffuse interstellar clouds (with masses 
estimated independently from observations of H I and CO) allow the opacity k of interstellar dust 
to be determined. The observed far-infrared and submm emission from diffuse clouds is consistent 
with f3 « 1.7 (Finkbeiner, Davis, & Schlegel 1999; Li & Draine 2001). 

The opacity spectral index /? can also be determined for dust in dense molecular clouds. The 
dust near embedded infrared sources in the Orion ridge has /3 « 1.9 between llOO/zm and 450//m 
(Goldsmith, Bergin & Lis 1997). The dust in and around 17 ultracompact H II regions studied by 
Hunter (1998) has (3 ^ 2.00 ± 0.25 between 350 and 450|Ltm. Priesen et al. (2005) find (3 = l.Qtoi 
for the dust in 3 hot molecular cores, with densities nu 10^ cm~^. The dust in both diffuse clouds 
and dark clouds thus appears to be consistent with an opacity spectral index /^ism 1.8 ± 0.2 in 
the submm. 

If the dust in protoplanetary disks were both optically thin and similar to interstellar dust, 
we would expect a = 2 + f^ism ~ 3.8 ib 0.2. However, Beckwith and Sargent (1991) found that 
protoplanetary disks usually have spectra with 2 < a < 3. The opacity spectral index /Sdisk was 
estimated by fitting simple disk models; the 24 objects studied had a median /3disk = 0.92, with 2/3 
of the sample falling in the interval /3(jisk = 0.9 it 0.7. 

Numerous subsequent studies have also found /3disk < Asm- Natta et al. (2004) find a median 
/^disk = 0.7 for a sample of 9 disks around Her big Ae stars observed with the VLA. Submm 
interferometry of a rotating disk in the star-forming region IRAS 18089-1732 indicates a — 2 3.2 
in the outer envelope, with small-scale structure having a — 2 ~ 0.5 (Beuther et al. 2004). Andrews 
& Williams (2005) obtained submm SEDs for 44 protoplanetary disks in Taurus-Auriga, and found 
that they can be fit using disk models with /9disk ^ 1- 

There are several possible explanations for the difference between Asm and /^disk^ 

1. Some of the emission may come from regions that are not optically thin, even at mm wave- 
lengths. This appears to be part of the explanation, but radiative transfer models also require 
material with /3disk ~ 1 in order to reproduce the observed spectra (e.g., Beckwith &; Sargent 
1991; Natta et al. 2004; Natta & Testi 2004; Andrews & Williams 2005). 

2. The chemical composition of the particles in the disk might be very different from the com- 
position of interstellar grains. 

3. Growth by coagulation in protoplanetary disks might lead to "fluffy" structures, with the 
difference between Asm and /3disk being due to differing grain geometry. 

4. As has been proposed by a number of authors (e.g., Beckwith & Sargent 1991; Miyake &: 
Nakagawa 1993; D'Alessio et al. 2001; Natta et al. 2004; Natta k Testi 2004) grain growth 
might have put an appreciable fraction of the solid mass into particles that are so large that 
they are not in the Rayleigh limit, even at mm wavelengths, leading to reduced /?. 
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Here we show that power-law size distributions dn/da oc for a < Omax, with 3 < < 4 and 
flmax ^ 3 A, will lead to opacity spectral index /3 ~ (p — 3)/5s, where (3s is the opacity spectral index 
of the solid material in the Rayleigh limit. Therefore, starting with material with the observed 
Pa ~ 1-7 of interstellar dust, grain growth with the likely power-law index p k, 3.5 will produce 
a mixture of sizes with /? 0.9 if Omax ^ 3A. Since grain growth is expected to be rapid in 
protoplanetary disks (e.g., Beckwith et al. 2000; DuUemond & Dominik 2005), this can naturally 
explain {3^1 observed for protoplanetary disks, without need to appeal to changes in composition 
or grain geometry, other than size. 

We briefly review what is expected regarding the optical constants of candidate materials at 
submm wavelengths in §2, and conclude that changes in composition are not likely to acccount 
for the low [3 values in disks. Particle size distributions are briefly reviewed in §3, to motivate 
the adopted power-law distribution. In §4, the opacity is estimated analytically for power-law size 
distributions, to find the expected behavior. In §5 we calculate the opacity numerically for several 
candidate materials. Although the dependence on size and wavelength can be complicated for 
conducting materials (in particular, amorphous carbon), we can reproduce the observed /3disk ~ 1-0 
for candidate materials with very different optical constants. Our results are discussed in §6 and 
summarized in §7. 

2. Dielectric Functions: What is Expected for /^^T 

Let /3a-(A) = d\a n/ dhxp in the small-particle limit, where n is the absorption cross section per 
unit mass of solid material. Simple models of insulators and simple models of conductors both have 
/3s = 2 at low frequencies (Draine 2004). How do real materials behave at submm wavelengths? 

If the low values of f3 observed in protoplanetary disks are to be produced by materials with 
small values of /J^, this would require materials with enhanced low-frequency absorption. In insula- 
tors, this would require optically-active low energy modes. Amorphous solids, with low-temperature 
specific heats revealing a spectrum of low-energy excitations (Pohl &; Salinger 1976), are candidates 
for enhanced low-frequency absorption. 

Laboratory measurements of the opacity of amorphous insulators are shown in Figure 1. In 
some cases (e.g., MgO-2Si02; Agladze et al. 1996) the submm opacity is consistent with a power 
law over the wavelength range studied. In other cases, broken power laws provide an acceptable 
fit. However, in all cases the power law index fig > 1- The smallest effective value is [3s ~ 1.2 for 
MgO-2Si02 between 150 and 300 GHz (Agladze et al. 1996), but for most samples, Ps ~ 2, with 
13s > 2 seen in some cases (e.g., Na20-3Si02 between 300 and 1200 GHz: Bosch 1978). While the 
number of materials studied is limited, it appears that small particles of amorphous insulators are 
likely to have (3s > 1.5 in the submm. 

For carbonaceous materials the situation is less clear. Opacities measured by Mennella et 
al. (1998) for materials "BE" (soot produced by burning CqB.q in air) and "AGAR" (amorphous 
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Fig. 1. — Laboratory measurements at various temperatures of opacities of amorphous Fe2Si04 (Mennella et al. 
1998), Na20-3Si02 (Bosch 1978), and MgO-2Si02 and MgO-Si02 (Agladze et al. 1996). Curve labelled WDOl is the 
opacity of the Weingartner & Draine (2001) dust mixture, for comparison. Dotted lines show power laws with (i = 1, 
1.5, and 2. 

carbon grains produced from a carbon arc in Ar gas) are shown in Figure 2. The opacities are very 
large (e.g., at 1mm, k « 25 — 30 cm^ g~^ for BE and AGAR at 24K, a factor 20-100 larger than for 
the insulators in Fig. 1 at T < 30 K). Samples BE and AGAR are characterized by small (5s'- at 
24K, BE and AGAR have (3s ~ 1-1 and 1.2, respectively. There is reason for concern that the very 
large submm opacities seen in the lab for materials BE and AGAR may be the result of contact 
between the particles in the sample (which could dramatically affect the opacity for materials with 
|e| ^ 1) in which case the measured opacities do not correspond to those that would be produced 
by isolated small spheres of the same material. 

Jager et al. (1999) have produced amorphous carbonaceous solids with varying levels of graphi- 
tization by pyrolysis of cellulose, and measured the optical constants. Pyrolysis at 600G in Ar gas 
produces a fine-grained material "cel600" with G:H:O::l:0. 33:0. 053 and density p = 1.67gcm~^. 
Although the material is substantially "graphitized" (in the sense that ~ 70% of the G atoms 
appear to be sp^ bonded, as in graphite), the material has 62 ^ as — > 0, implying zero d.c. 
conductivity. Jager et al. provide the refractive index for 0.2/im < A < 450/im. We extrapolate 
to lower frequencies assuming ei(A > 450/im) ~ 4.03 and 62 (A > 450//m) = 0.0378(450/im/A)^'^, 
giving (5s = 1-3 for A > 450/Lim. The opacity calculated for small spheres of this material is shown 
in Figure 2. As noted by Jager et al., small spheres have (5s ~ 1-5 for A > lOO/zm. 
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Fig. 2. — Same as Figure 1 but for carbonaceous materials BE and AGAR measured at T = 24K and T — 295 K 
(Mennella et al. 1998), and calculated for small spheres with optical constants measured for cellulose pyrolyzed at 
600C and 800C (Jager et al. 1999) - see text. 

Pyrolysis at 800C produces a material "celSOO" with C:H:O::l:0. 14:0.030, with p = 1.84 gcm-^ 
(Jager et al. 1999). The estimated sp'^ fraction is ~ 75%, and the material becomes conducting. 
Jager et al. measured the optical constants for 0.178/um < A < 562/xm; we extrapolate to lower 
frequencies with a simple free-electron model with conductivity a = 3.6 x 10^^ s~^, and damping 
time T = 1.5 x lO^^^s. The opacity calculated for small spheres of cel800 is shown in Figure 2. 
At 1 mm, this calculated opacity is more than an order of magnitude below what was reported for 
AGAR and BE. As noted by Jager et al., Ps = 1-95 for A > 100/im. 

With |e| ^ 1 in the submm, k for cel800 is strongly shape-dependent. Because of its simple 
analytic properties, the "Continuous Distribution of Ellipsoids" (CDE) is sometimes used to allow 
for nonspherical shapes. Jager et al. show that for the CDE, celSOO would have (3s ~ 0.7, but 
caution that the CDE serves only as an illustrative example of shape effects. The CDE assumes 
that a small fraction of the ellipsoids have very extreme shapes, and needle-like particles make a 
large contribution to n when |e| ^> 1 (Min, Hovenier, &: de Koter 2003). The shape distribution 
of interstellar particles does not seem likely to resemble the CDE, at least as regards these ex- 
treme shapes. We consider that spheres (or ellipsoids with modest axial ratios) provide a better 
approximation to interstellar particles than the CDE. 

Any material capable of explaining /3 < 1 on the basis of composition rather than size would 
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have to be unlike the insulators in Figure 1, or the pyrolyzed cellulose in Figure 2. The opacities 
reported by Menella et al. for carbonaceous materials BE and AGAR are large and have /3 ~ 1, 
but it is not clear that the measured opacity would apply to free particles and, if it did, such 
material could not be a significant component of interstellar dust. Although one cannot exclude 
the possibility that protoplanetary disks might include material with both large mm/submm opacity 
and /3s 1, alternate explanations should be considered. 



3. Size Distribution 

Consider particles with a power-law size distribution, 

dn . f a ^ ^ 

— = A \ amin <a< Omax , (3) 

da VOmax/ 

where A is a constant; the size distribution presumably results from competition between growth 
(e.g., coagulation) and destruction (e.g., shattering). 

Experimental studies of fragmentation find power-law size distributions for the fragments, with 
p varying between ~ 1.9 for low- velocity collisions, topK. A for catastrophic impacts (Davis & Ryan 
1990). However, the power-law index for fragments from shattering events integrated over a range 
of target masses need not be the same as the power-law index for fragments of a single target. 

Dohnanyi (1969) showed that, for certain assumptions, coagulation and coUisional fragmenta- 
tion would lead to a steady-state size distribution with p = 3.5. Tanaka et al. (1996) argue that 

the p = 3.5 power law is in fact a very general result that depends only on the assumption that the 
fragmentation process is self-similar, and that the collision rate varies as a? . Self-similarity may 
not apply to larger bodies, where self-gravity becomes important, or to the smallest sizes, where 
surface free energy considerations will limit the numbers of very small fragments, but self-similarity 
may apply, at least approximately, to the cratering and shattering events occuring in the interstellar 
medium and in protoplanetary disks. 

Extinction and scattering of starlight by interstellar dust can be reproduced by a mixture of 
graphite and silicate grains with p k. 3.5 and Omax ~ 0.25/xm (Mathis, Rumpl, k. Nordsieck 1977; 
Draine & Lee 1984). Although more recent work (e.g., Kim & Martin 1995; Weingartner & Draine 
2001) find deviations from a simple power law with cutoffs, the p = 3.5 power law gives a good 
accounting for the overall size distribution spanning many orders of magnitude in mass. While 
the origin of the interstellar grain size distribution remains uncertain, it is intriguingly close to the 
predictions of self-similar fragmentation. 

The distribution of D > 30 km asteroids has p 3.4 (Dohnanyi 1969). Based on observations 
by the Sloan Digital Sky Survey (Ivezic et al. 2001), the asteroid size distribution appears to be 



^If the collision rate varies as a', then p = 2.5 -|- q/2 (Tanaka et al. 1996). 
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consistent with p 3.25 for diameters D between 5 and 300 km (Bottke et al. 2005). 

If the size distribution of particles in a protoplanetary disk is determined by a competition 
between coagulation and fragmentation, it is plausible that the size distribution might be close to 
a power law with p ~ 3.5, with an upper cutoff amax that will depend on how far coagulation has 
proceeded. Note that for p ~ 3.5, most of the mass is at the large size end of the distribution, with 
the numbers of large grains nearly independent of the lower cutoff Omiin provided Omin ^ cimax- 



4. Analytic Approximation for the Opacity 

For a size distribution dn/da, the opacity (absorption cross section per unit mass of solid 
materials) at frequency v is given by 

^ / da{dn/da)Cahs{a,i^) ... 

~ fda{dn/da){ATr/3)pa^ ' ^' 

where p is the solid density, and a = (3m/47rp)^/'^ is the effective radius of a particle of mass m. 

Consider some frequency v, wavelength X = c/v, and suppose that the absorption cross section 
for particles of size a can be approximated by 



{A7r/3)a^a,y for a < adi') 
'Kc? for a > ac(z^) 



where 



IStt £2 , _ 3 A (ei + 2)2 + ej 

Oi/ = t:^ o a^nd Uc = - — = — (6) 

A (ei + 2)2 + ei 247r £2 ^ ^ 

are defined so that eq. (5) is exact for spheres in the small particle limit (a/A — 0) (see, e.g., 
Draine &; Lee 1984). Here ei(A) and £2 (A) arc the real and imaginary parts of the complex dielectric 
function e = ei + ie2- Eq. (5) is only approximate for a > X/2tt, but the limiting behavior in the 
geometric optics limit a ^ cx) is sensible, corresponding to a low-albedo object. Because we will be 
interested in size distributions where the surface area is dominated by small particles, our analytic 
estimate of the opacity integral (4) is relatively insensitive to errors in the albedo of the particles 
with a S> A/27r. For a ^ A/27r, Cabs can exceed eq. (5) by as much as a factor ~ 10 (depending on 
e), but for the moment we use eq. (5) to estimate the opacity for the size distribution (3): 

k(i/) = ^ if Omax < Oc (7) 

and 

4—p / \ 3—p 

1-{P-S)(~] -i4-p) 



P (p-3) 



Or. 



(8) 



for Omin < Qc < Omax, p > 3, p ^ A, and amin < flmax- If 3 < p < 4, then 

3 1 

4/9aniax {P - 3) 



'^M ~ 1 7 ^ o " Omin < flc < Omax • (9) 
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If the small-particle absorption cross section per volume varies as a power-law in frequency, 

a^ = aQ\ — \ , (10) 

then 



'^('^) = 7 ^ 1 ~ if "min < ac <C amax ■ (H) 

Therefore, if Omin <C Oc <C amax, the size distribution will have 

P^ip-3)Ps . (12) 

For p = 3.5 and Ps ~ 1-7, eq. (12) predicts /? ~ 0.85: material with the FIR-submm opacity 
characteristic of diffuse interstellar dust (/3s ^ 1.7) can reproduce the observed frequency depen- 
dence of the opacity of material in circumstellar disks (/3disk ~ 1) if cimax is large enough so that 
Omax S> ac{X) over the frequency range of interest. 

The derivation of eq. (11-12) relied on the approximations (5), which are known to not be 
quantitatively accurate. We now proceed to direct calculation of /^{v) for several candidate materials 
to determine the actual frequency dependence of (3, and to see what values of k result when amax 
is large enough to satisfy the condition (9). 



5. Examples 

The precise behavior of k(z/) depends on the complex dielectric function e(z/). Not knowing the 
precise composition of interstellar and protostellar grain material, we consider three, quite different, 
candidate materials. We will see that for the different materials considered here, a power-law size 
distribution with p 3.5 results in mm/submm opacities with (3 ^ 0.9 provided Omax ^ 3 mm. 

The opacity k{i') (eq. 4) is evaluated with Cabs(a) ^) calculated for spheres, using Mie theory for 
a < Qx = 2x 10^A/27r|m|, and geometric optics for a > Qx, where m = ^/e is the complex refractive 
index. We employ a double precision version of the Mie theory implementation of Wiscombe (1980, 
1996). 



5.1. Silicate Particles 

We calculate the opacity of spherical silicate particles using the dielectric function estimated 
for interstellar silicate material (Draine 2003). Figure 3 shows k(i/) for size distributions with 
p = 3.5, amin = 3.5 A, and various values of Omax- The curve labelled 0.25//m is the "small 
particle" opacity estimated for interstellar silicate grains. At A = 1mm, increasing amax does 
not noticeably change the opacity until amax reaches 100/xm, at which point the optics of small 
particles with 27ra/A ^ 1 results in an increase in the opacity. This effect leads to opacities for the 
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overall size distribution that exceed the small-particle opacity at 1 mm until amax reaches 10 cm. 
For Oinax ^ 1mm, K(lmm) oc l/^fflmax, consistent with eq. (11). The upper panel in Fig. 3 shows 
the opacity spectral index /3(A) for the different choices of Omax- For this material, the condition 
ac <C Cmax can be replaced by 3A < Cmax- We see that /?(lmm) < 1 for Omax ^ 3 mm. 
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Fig. 3. — Opacity of amorphous silicate spheres with size distribution dn/da oc a~'^'^ for 3.5 A < a < Omax- Curves 
are labelled by flmax; curve for 0.25/im is the small-particle limit, and numerically close to the WDOl dust model (not 
shown). Upper panel shows /? = din n/dlnv for selected flmax- /3(lmm) < 1 is found for Omax 3mm. 



5.2. Amorphous Carbonaceous Material 

The dust model of Weingartner & Draine (2001; hereafter WDOl) includes carbonaceous mate- 
rial with the anisotropic dielectric function of graphite, with absorption and scattering by randomly- 
oriented spherical grains estimated using the "1/3-2/3 approximation" (Draine & Malhotra 1993). 
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Fig. 4. — Same as Fig. 3, but using the optical constants of cellulose pyrolyzed at 600C (Jager et al. 1999). Also 
shown is the opacity for the WDOl dust model. For Omax < 100/im, K(lmm) is very similar to the WDOl opacity. 
For amax }^ 3 mm, the opacity has P{1 mm) < 1. 



Graphite is a conductor, with appreciable conductivity in the basal plane, and nonzero electri- 
cal conductivity even parallel to the c axis (i.e., normal to the basal plane). For large particles 
at submm wavelengths, the conductivity plays a major role in determining the absorption cross 
section. Because large particles formed by coagulation will not have the properties of crystalline 
graphite, we consider instead the amorphous carbonaceous solids celGOO and celSOO studied by Jager 
et al. (1999) (see §2). 

Figure 4 shows the FlR-submm opacity for size distributions of celGOO spheres. In the small- 
particle limit, the opacity of cel600 near A = 1 mm is similar to or less than the opacity of the WDOl 
grain model, and at 100;um the opacity is less than the WDOl model by a factor ~ 3. The WDOl 
model successfully reproduces the observed emission from the ISM (Li & Draine 2001) and other 
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spiral galaxies (Regan et al. 2004). Candidate interstellar grain materials that do not substantially 
exceed the WDOl opacity could, in principle, be abundant in the ISM. Thus, the mm/submm 
properties of interstellar carbonaceous material could perhaps resemble celGOO. The mm/submm 
opacity of the size distributions of celGOO is insensitive to Cmax until Omax reaches ~ lOO/im. Size 
distributions of celGOO with p = 3.5 have /3(1 mm) < 1 for Omax ^ 1 cm. 




Fig. 5. — Same as Fig. 4, but for cellulose pyrolyzed at 800C (Jager et al. 1999), resulting in a material with 
appreciable d.c. conductivity. In this case, magnetic dipole absorption (see text) causes k(A — 1 mm) to increase as 
fflmax is increased beyond Ifim. These particles have /3(lmm) < 1.2 for flmax ^ 200^m. 

Spheres composed of celSOO, on the other hand, gives FIR and submm absorption that is 
much stronger than in the WDOl dust model, even when the material is in the form of spheres - 
any other shape would provide even stronger absorption. For 1 mm ^ A < 100/im, in the small- 
particle limit K exceeds the WDOl opacity by factors of ~ 4; hence it does not seem likely that 
the interstellar grain mix could include very much such material. Nevertheless, because we cannot 
rule out the possiblity that protoplanetary disk material might have optical properties resembling 
celSOO, we consider size distributions of spheres composed of celSOO as an example of a conducting 
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Fig. 6. — /3(j^) for size distributions dn/da oc , for 3 values of p, and the materials of Figs. 3-5. Dotted lines 
are (p — 3)Ps, which should be approximately equal to /3(A) (see eq. 12), for p — 3.2, 3.5, 3.8 (bottom to top in each 
panel) showing approximate agreement with the actual /3 values. 



material. In the small-particle limit, celSOO has /3<j ~ 2, but size distributions with Cmax ^ 3 mm 
give /3(lmm) ~ 0.9, showing the generality of the result (12). 

To see the sensitivity of /5 to the power-law index p, Fig. 6 shows /? for p = 3.2, 3.5, and 3.8. 
The opacity spectral index (3(1 mm) < 1 requires p < 3.6 for astrosilicate or celGOO, or p < 3.8 for 
cel800. 

The relation f3 (p — 3)Ps was derived using eq. (5), which is not an accurate approximation 
for Cabs when a > Xjl-n. Also shown in Fig. 6 (dotted lines) are (p — 3)/9s(z^) for p = 3.2, 3.5, 
and 3.8. We see that while (p — 3)/?^ is often close to the actual there are regimes where /3 

differs appreciably from this estimate (e.g., astrosilicate with p = 3.8 at > 10^ GHz; celGOO with 
p = 3.2 for u < 10^ GHz, and cel800 with p = 3.8 for u < 500 GHz). Nevertheless eq. (12) provides 
a reasonable estimate. 
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6. 



Discussion 



6.1. 



Ices? 



In diffuse clouds, H2O and other ices are not detected: ultraviolet radiation prevents the ac- 
cumulation of H2O ice coatings on grains in diffuse clouds. In dark clouds, however, substantial 
amounts of solids containing H2O, CO. CH,30H, H2CO, NH3, CH4, and other compounds are ob- 
served to be present, presumably as "mantles" coating "cores" (silicate and carbonaceous material) 
inherited from diffuse clouds. 

What will be the role of these ices at submm wavelengths? Using the dielectric function of 
crystalline H2O ice I, Aannestad (1975) found that the presence of ice mantles had only minor 
effects on the absorption at A > 200/xm, because the vibrational and hindered-rotation modes of 
the H2O all lie at higher frequencies. This is confirmed in recent work (Pontoppidan et al. 2005). 

Preibisch et al. (1993) considered silicate grains coated with "dirty ice" consisting of a mixture 
of H2O and NH3 with small amorphous carbon inclusions. The optical properties of this mixture 
were estimated using effective medium theory. In the submm, the absorption by this "dirty ice" 
material was almost entirely due to the amorphous carbon inclusions. When present as inclusions 
within an ice matrix, the amorphous carbon material is more absorptive than when present as a bare 
particle. It should be kept in mind, however, that this conclusion is based on use of approximate 
effective medium theories to estimate the optical properties of the ice with inclusions. 

Therefore, although the absorption cross section for a silicate or carbonaceous grain will change 
if an ice mantle develops on the grain, at submm frequencies it does not appear that this will 
substantially alter the grain absorption cross section unless the ice itself is strongly absorbing, and 

strong absorption does not appear likely unless the "ice" is contaminated by amorphous carbon 
particles - in which case the absorption is due to the amorphous carbon, not the ice itself. 



Beckwith &: Sargent (1991) suggested that grain growth might explain the values of /3disk ~ 
1 inferred from their observations. Miyaka &: Nakagawa (1993) calculated the opacity for size 
distributions extending to Omax ^ 1 cm, and concluded that grain growth could indeed account for 
the observed values of /3disk- Miyaka &; Nakagawa considered grains composed of mixtures of silicate 
and H2O ice, either as compact solids or in fluffy structures. The present paper shows that the 
result of Miyaka &; Nakagawa is quite general, and applies to both conducting and nonconducting 
materials. 

D'Alessio et al. (2001) calculated the opacity for power-law size distributions with p = 2.5 and 
p = 3.5, for a mixture of separate spheres of silicate, ice, and organic material, but for the silicate 
material they assumed constant €2 for A > 800//m, resulting in /?s = 1 for the silicate spheres at 
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A > 800|Um. Therefore even in the small particle limit, their model had /?s = 1 for A > 800/im. 
This frequency dependence is not expected - see Fig. 1 where amorphous insulators have /?s 2 at 
mm wavelengths. 

Natta ct al. (2004) calculated k and (3 for power-law size distributions with various exponents 
p for mixtures of compact spheres of olivine, organic material, and water ice, and also for fluffy 
grains containing mixtures of the above, with 50% vacuum. They used dielectric functions such 
that /3s = 1.27 for the mixture of separate, compact spheres (Natta k, Testi 2004). They showed 
that by increasing Omax to Omax ^ 3 mm for the compact grains, or Omax ^ 10 cm for the fluffy 
grains, size distributions with p 3.5 would have /3(1 mm) 1. 

The results of the present paper are in general agreement with these previous studies. We also 
point out that when Omax ^ A, the approximate estimate /? ~ (p — 3)/3s applies to both insulating 
and conducting particles (for 3 < p < 4). For the likely power-law index p ~ 3.5, this gives 
(3 ~ 0.5/3s, so that materials with (ig = 2, for example, are expected to have /? ~ 1 when Cmax ^ A. 
Direct calculation of n for size distributions of three candidate materials - silicate, nonconducting 
carbonaceous material (cel600) , and conducting carbonaceous material (celSOO) confirm that power- 
law size distributions with p 3.5 will have /? < 1 for Omax ^ 3A. 

It is sometimes suggested that the low values of /?disk are due to "fluffiness" of grains grown by 
agglomeration (Beckwith & Sargent 1992). Fluffiness could subtantially increase the mm/submm 
opacity if the grain material is electrically conducting (Wright 1987), but for nonconducting mate- 
rials with dielectric functions that are not very large at the wavelengths of interest, grain geometry 
will have only modest effects on the grain opacity. Kriigel & Sicbcnmorgen (1994) considered ex- 
tinction and absorption by "fluffy" grains composed of silicate, amorphous carbon, and ice. The 
absorptivity of the ice in the far infrared was treated as a free parameter, and chosen to be large 
enough for the ice absorption to dominate the submm absorption. With their adopted optical 
constants, the particles had /3s 1.8 for A > 100/xm. They found that power-law size distributions 
(3) with p = 3.5 and amax = 3 mm had (3 ^ 0.5 for lOO/xm < A < 2 mm. 

Stognienko et al. (1995) discussed the optical properties of particles produced by coagulation 
of smaller sub-particles, and showed that the opacity of such clusters depends strongly on the 
optical properties of the sub-particles, and on the geometry of the clustering. For amorphous 
carbon sub-particles, the particle geometry can lead to strong enhancements in absorption for some 
cluster topologies. However, if the cluster size a <C A, the clusters appear to have /3 « 2 at submm 
frequencies. Henning & Stognienko (1996) calculated the opacity for dust aggregates; their standard 
model has /3 Ri 2 for the aggregates. 

From the above studies we conclude that, while grain structure unquestionably plays a role 
in determining scattering and absorption cross sections, the observed small values of /3disk are 
naturally explained simply as the result of increased grain size, independent of whether the grains 
are compact or conglomerate. Grain growth to sizes > 3 mm can naturally account for the observed 
/3disk(lmm) Ri 1. 
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7. Summary 

The principal conclusion of this paper is the following: if the solid particles in protostellar disks 
have size distributions dn/da oc for a < Omax, with p « 3.5 (as observed for interstellar dust, 
and as seems likely for particles in protoplanetary disks) and Omax ^ 3 mm, the resulting opacities 
have /9(lmm) « (p — S)/?^ « 1, where is the opacity spectral index in the small particle limit. 
Therefore if the particles in protoplanetary disks are composed of the same material as interstellar 
grains (with (3s ~ 1.8 it 0.2 in the submm) and the size distribution has p ~ 3.5 and extends to 
Omax ^ 3mm, the protoplanetary disk material will have /3(lmm) < 1. This result is not expected 
to be altered by the addition of of dielectric ice coatings on the grains. 

To demonstrate the generality of this conclusion, three quite different materials were consid- 
ered: two amorphous insulators (astrosilicate and cellulose pyrolyzed at 600C) and a conductor 
(cellulose pyrolyzed at 800C). While the relation /? ~ (p — 3)/3s is only an approximation, we find 
that all 3 materials give /3(1 mm) < 1 for p 3.5 and Omax ^ 3 mm. 

The observed submm emission from protostellar disks therefore is explained naturally by the 
particle growth that is expected to occur in these disks. 

This research was supported in part by NSF grant AST-0406883. I am grateful to A. Natta 
for helpful discussions, and to R.H. Lupton for availability of the SM software package. 
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